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Abstract 

We consider a model of DBI inflation where two stacks of mobile branes are moving ultra relativistically 
in a warped throat. The stack closer to the tip of the throat is annihilated with the background anti-branes 
while inflation proceeds by the second stack. The effects of branes annihilation and particles creation during 
DBI inflation on the curvature perturbations power spectrum and the scalar spectral index are studied. 
We show that for super-horizon scales, modes which are outside the sound horizon at the time of branes 
collision, the spectral index has a shift to blue spectrum compared to the standard DBI inflation. For small 
scales the power spectrum approaches to its background DBI inflation value with the decaying superimposed 
oscillatory modulations. 
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INTRODUCTION 



Inflation proved to be very successful as a theory of early universe and structure formation {l| 
which is strongly supported by cosmic observations (sj. The simplest models of inflation are based 
on a scalar field minimally coupled to gravity with a potential flat enough to support an extended 
period of inflation. Despite its observational successes there is no deep theoretical understanding of 
inflation and the nature of the inflation field. There have been many attempts during past decade 
to embed inflation within the context of string theory for a review see [3-^] and references therein. 



Brano inflation is an interesting realization of inflation Iron) string Ihoory [9l-l^|. In ils original 
form, it contained a pair of D3 and anti D3 branes moving in the Calabi-Yau (CY) compactification. 
The inflaton field is the radial distance between them so in this sense the inflaton field has a 
geometric interpretation in string theory. Inflation ends when the distance between the brane 
and the anti-brane reaches the string length scale where a tachyon develops in the open strings 
spectrum stretched between the pair. Inflation ends soon after tachyon formation and the energy 
stored in branes tensions are released into closed string modes However, it was soon realized 
that the potential between the pair of brane and anti-brane is too steep to allow a long enough 
period of slow-roll inflation. To flatten the potential, it was suggested to put the pair of brane and 
anti-brane inside a warped throat 



as in Randall-Sundrum scenario 



lc -41511 . where the potential between D3 and D3 is warped down 



UrkZl]. 



As a novel feature of brane inflation the question of branes annihilation and particles creation 
during inflation was studied in [^j]. In this picture two stacks of branes, located at different places 
inside the throat, are moving in the slow-roll limit towards the bottom of the throat where anti- 
branes are located. The stack closer to the tip will be annihilated during inflation transferring its 
energy into closed string modes. The second stage of inflation is driven by the remaining stack of 
branes until it is annihilated by the remaining anti-branes at the bottom of the throat. The process 
of fields annihilation and particle creations 



31 



371 ] can have interesting observational 



n 

As in [28J], the stack closer to 



consequences on CMB. Here we generalize the results in [281 ] to_the case where the mobile stacks 
of branes are moving ultra relativistically as in DBI inflation 
the tip is annihilated by the background anti-branes resulting in particles creation during inflation 
while the second stage of inflation is driven by the remaining stack. 

The rest of the paper is organized as follows. In section [II] we present our set up and background 
solutions for DBI inflation. In section IIIII we study the curvature perturbations and present our 
matching conditions. In section IIVI we obtain the power spectrum of curvature perturbations and 
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FIG. 1: A schematic view of the set up. There are two stacks containing p\ and P2 coincident branes inside 
the throat. The stacks are moving relativistically towards the bottom of the throat where there are p\ +P2 
anti- branes. The figure is borrowed from [3] 

the spectral index. Discussions and conclusions are presented in section [V] 

II. THE MODEL AND BACKGROUND EQUATIONS 

In this section we present our set up. As explained above we consider two stacks containing p± 
and P2 coincident branes. The positions of the stacks are represented by functions r\(t) and ^(t) 
inside the throat. There are p\ +P2 anti-branes located at the bottom of the throat, ro- We assume 
that the first stack is closer to the tip of the throat so it is annihilated during inflation transferring 
its energy into closed strings modes. The second stage of inflation is driven by the remaining stack 
of P2 branes. For a schematic view, see Fig. Q]and Fig. [2j 

As usual the metric of the throat is 



where r is the radial coordinate, dQ,^ represents the internal five-dimensional azimuthal directions 




(1) 



which we may take to be an S 5 and h(r) is the warp factor 




(2) 
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FIG. 2: A schematic view of the annihilation event at t c . The stack closer to the tip of the throat is 
annihilated during inflation transferring its energy into radiation. The second stage of inflation is driven by 
the remaining stack containing p2 branes. Inflation ends when this stack is annihilated by the remaining p2 
anti-branes at the bottom of the throat. The figure is borrowed from 281 ] . 



Here L is the AdS length scale of the throat which is created by N coincident background branes 
located at r = (for a review see 39^ ) . 

L 4 = Ang s Na' 2 , (3) 

where g s is the perturbative string coupling and a' = I 2 where l s is the string length scale. The 
mobile stacks of p\ and p 2 branes are located at positions r\ and r 2 as probe branes. In order for 
our probe branes approximation to be valid, we need pi,P2 <C N. In addition there are p\ +P2 D3 
branes at the bottom of the throat ro- Once the distance between the stack of p\ branes and these 
anti-branes reaches at the order of l s the open strings stretched between them become tachyonic, 
indicating an instability in the system. Shortly after tachyon formation, the stack of p\ branes 
is annihilated by p\ anti-branes at the bottom of the throat. Consequently, the energy stored 
in pi branes and anti-branes tension, which is 2p\T^h{ro)~ 2 , is transferred into closed strings 
modes. For simplicity, we assume that these closed strings modes are in the form of massless 
particles (radiation). Of course a combinations of massive and massless particles would be created. 
However, once these particles are created they will be diluted quickly by the background inflation 
so for our practical purposes, there are not much differences between massive and massless particle 
creations. 

The goal of this work is to examine the effects of branes annihilation and particles creation on 
the remaining field which drives the second stage of inflation. In 28] this was studied for the case 



where the stacks are moving slowly. Here we genera 
the stacks are moving ultra relativistically 



38 



40H44]. 



ize those analysis to our case at hand where 



The Dirac-Born-Infeld (DBI) action for the stacks of p\ and pz branes moving relativistically in 
the background of Eq. dU is 

S = E^ T 3 / V=gd 4 x h~\n) (l - y/l - h(rj) rjj - Vi(n) , (4) 



where T-* is the D3-brane tension with T^L 4 



N 



We also have added the unknown contributions 



'.3 is tne uo-orane tension wren , . 

V\(r\) and V-zir-z) from the back-reactions of the background fluxes, branes and Kahler modulus 
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221. As usual, we shall continue from the 



38, 



stabilization on the dynamics of mobile branes 

phenomenological point of views and parametrize these potentials to be suitable for inflation 

To proceed further, we simplify the DBI action into the following form which is more suitable 
for our analysis 

2 



1=1 J 



xv-g 



(5) 



where 0/ = y/T^piri, fi=ty and A/ = T^L^pj = -^pf. This form of multiple field DBI inflation 



has been studied extensively in the past 



47J 



Now we are ready to promote the system into a cosmological set up. The background four- 
dimensional FRW metric is 



ds z = -di A + a{tfdx 2 , 



(6) 



where a(t) is the scale factor. The Friedmann equation and the energy conservation equations are 

(7) 



3# 2 = T72 > p + 3H(p + P) = 0. 



Ml 

where H = a/a is the Hubble expansion rate. The energy density p and the pressure p are given 
by 



p=T / iK\n-i)+vi] 



P = E[/T 1 (l-77 1 )-Vr] 



(8) 



I 1 
Here the Lorentz factor for each field, 7/, is defined by 

1 



1 - flfl 



(9) 



Finally, the background Klein-Gordon equation for each field is 

h + aff 7 rVi + ~fi ft ~ | + 77 3 (v/ + 1) = o . (io) 

So far we have not specified the form of the potentials, Vj ((f) j)- As explained above, these 
potentials arise from the back-reactions of the fluxes, branes and Kahler modulus stabilization to 
the mobile branes. It is a non-trivial question in string theory as how one can calculate these 
back-reactions 



23, 



24 



26[ concretely. However, to be specific, we consider the phenomenological 



+ -75 = °, (11) 
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approach where the potentials have the form Vj = -^4>f- We shall also keep the masses mj 
undetermined except that their magnitudes should be such that they give the desired number of 
e-foldings and fit the CMB observations such as the COBE normalization. 

We are interested in the limit where the branes are moving ultra relativistically inside the throat 
so one can not expand the square root in Eq. ([5]) perturbatively in terms of <f). This corresponds 
to the case where 7 > 1 and the stacks of branes are approaching the speed limit 0| — fj . In 
this limit Eq. (fTUl) reduces to 

which can be verified to have the speed limit solution <pj(t) ~ ^y 1 - One can plug this solution 
back into Eq. (|10p to find the next leading correction and 

To get this solution, we have neglected the variation of H during inflation to leading order. 

To obtain the solution for the scale factor a(t) we require that the energy density is dominated 
by the potential energy so the Friedmann equation reduces to a/a oc 1/t. Below we find the 
background solutions for both inflationary stages, before the collision and after the collision. 

1. Before Collision 

Before collision both stacks contribute to the energy density and the Friedmann equation can 
easily be solved to give 

where the dimensionless number A_ is defined by 



Working with the number of e-foldings N(t) as the clock, dN = Hdt, one obtains 

JV(t) = A_ln-^, (15) 

where ti n is the time at the start of inflation when N = 0. Define N = N c the time when the first 
branes collision take place and the stack of p\ branes is annihilated by p\ background anti-branes. 
Similarly, define (pi = <j)i c as the corresponding fields values at the time of collision. Combining 
Eqs. (USD and ® yield 

N c 

4>ic ^ 4>n n e A - , (16) 

where (fin n are the corresponding fields values at the start of inflation. 

One can find an approximate formula for the time of branes collision as follows. As mentioned 
before, the time of branes annihilation is when the physical distance between the stack of p\ 
branes and anti-branes located at the bottom of the throat ro reaches the string length scale l s . 
Starting with the metric (H|) the physical distance between the first stack and the anti-branes, A£, 
is calculated to be 

/"T'l 

A£= h{rf/ A dr = L In — . (17) 
Setting A£ = l s and using Eq. (fT6j) the onset of branes collision, iV = N c , is obtained to be 

N r ~ A_ I In hm ls 



4>iA L 

^ A_(ln^-(4ng s N)-^ , (18) 

where to get the final relation Eq. ([3]) has been used. Here we defined (friA = VpiT^ro as the values 
of the corresponding fields at the position of the anti-branes. 

Correspondingly, the position of 4>2 at the time of collision is obtained to be 

VX~2 



2c ^ —— ^ fc»e A - . (19) 



Also it is instructive to calculate 77. Using Eqs. (|12|) and (fT3|) in the expression of 77- one finds 



m? T t 2 



77 = -at- (20) 



This indicates that as inflation proceeds, 7/ increases like t 



2 
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2. After Collision 



After the first stack is annihilated, the second stage of inflation is driven by the remaining stack 
of p2 branes and the scale factor is given by | = ^j- where now 

Defining the total number of e- foldings by Nt, with ~ 60 to solve the flatness and the horizon 
problem, one has 

N T - N c ~ A+ln^ , (22) 

where tt is the time of end of inflation when the second stack is annihilated by the remaining p2 
anti-branes. As before, this happens when the physical distance between the second stack and the 
anti-branes reaches the string scale l s which can be used to determine the value of 02 at the end 
of inflation 02/ — 4>2A exp(/ s /L). Similarly, one can show 

Nt ^ N c + A + ln — 

02/ 



- i i - ^ I a; + a + 1h 



(4TTg s NY 



(23) 



_4>2A 

One should arrange the background parameters such as g s N, A± and 0/ m such that large enough 
Nt can be obtained to solve the flatness and the horizon problem. 
Correspondingly, the Lorentz factor after the collision is 

+ [2m 2 M P t 2 

72 "V3^r- (24) 

One can find a measure of the magnitudes of A± as follows. For DBI inflation to match the 
COBE normalization, one requires a large background charges such that g s N 3> 1 {32I, 



61(] so 



one can safely neglect this term in Eqs. (J23J) and f 1 1 8 1) . Assuming that there is not exponential 
hierarchy between 02m an d 4>2A ° ne concludes that A + < Nt- However, by keeping the ratio 
exponentially large in the light of one can lower A + by a factor of ten or so. In our analysis 
below, we shall take take A + ~ A_ 3> 1. 
As shown in 



significant amount of non-Gaussianities can be created in DBI inflation where 
the branes are moving ultra-relativistically. The non-Gaussianity parameter /jvl is related to the 
Lorentz factor via /nl — — O.37 2 . To satisfy the WMAP constraints on f^L one requires that 
7 < 31 



41 



43j. 
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III. PERTURBATIONS 



Having studied our inflationary background in some details, now we consider the perturbations 
in this background. As discussed in 28[ a complete treatment of perturbations in our set up with 



branes annihilation and particles creation is a non-trivial task. The process of field annihilation 
and particles creation are determined by the dynamics of open string tachyon condensation. In the 
presence of background branes and fluxes this is a non-trivial phenomena 48|, [49( . We shall instead 
take the phenomenological approach and impose some physically reasonable approximations to 
bypass these difficulties. 

Before the collision we have two scalar fields 4>\ and <p2 corresponding to two stacks of branes 
moving relativistically. After the collision, the energy associated with field is transferred into 
closed strings modes, which we approximate them as massless particles. The process of field 
annihilation and closed strings formation takes some time scale, At. In our set up where the 
collision happens at the bottom of the throat At is given by the inverse of the warped string mass 
scale, that is At ~ /i(ro) 1//4 mJ 1 where m s = Ij 1 is the string theory mass scale. Note that the 



161 ] . In our effective four-dimensional 



factor /i(ro) 1//4 in At is in the light of Randall-Sundrum idea 
approach where H <C m s , we can take the process of field annihilation and particles creation as 
instantaneous and set At ~ for practical purposes. This approximations is certainly true for 
large scale perturbations. 

To be slightly more specific, we can borrow the results of 50-53]. As mentioned before, after 



brane and anti-brane annihilation their energy is transferred into closed strings modes. Subse- 
quently, these modes cascade down to graviton zero mode and the Kaluza-Klein (KK) excitations. 
This happens in a time scale controlled by the inverse of the warped sting mass scale, much shorter 
than H~ x . Subsequently, the massive KK modes decay into lighter modes. Because of the warp 
factor, the wave function of these modes are highly peaked at the bottom of the throat where the 
annihilation takes place. Their effects on the remaining stack of p2 branes at position r<i may be 
negligible in the limit where the stack of p2 branes is far away from the tip of the throat. 

Even in the instantaneous brane annihilation approximation it is not clear within string theory 
set up as how to turn off field (/>ijmd turn on radiation. To bypass this shortcoming we follow the 
practical approach employed in 



281 ]. We assume that field 0i is subdominant in the background 
energy density, corresponding to V± <C V<i- After §\ is annihilated, we assume that not only its 
energy is transferred into radiation energy density, p r , but also the perturbations carried by 54>i are 
transferred into perturbations in radiation, 5p r . We also assume that the perturbations 5p r , before 
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significantly diluted by the background expansion, do not interfere with 5 fa after collision. This 
can be considered as a higher order back-reaction, perhaps controlled by the ratio P1/P2, which we 
take to be small. Also this back-reaction may be negligible as long as the separation between two 
stacks of branes are large enough inside the throat. Within these approximations we do not need 
to follow the perturbations carried by 5<j)\ before the collision and the perturbations carried by 
5p r after the collision. In these approximations we shall follow only the perturbations 5(j>2 before 
and after the collision. This may not be such a bad approximation because physical quantities, 
such as curvature perturbations, are calculated at the end of inflation when only field fa and its 
perturbations are present. Since we do not follow 8<j)i and 8p r , our treatments of perturbations 
eventually reduce to those of single field model where only 5(j)2 is counted before and after the 
branes collision. However, the effect of (pi and radiation is kept in the background so that is how 
S<f>2 feels their presences. 



A. Perturbations Equations 

With these approximations we start the analysis of perturbations in our set up. The metric 
perturbations in the conformal Newtonian gauge where there is no anisotropy in stress energy 
tensor is 

ds 2 = -(1 + 2&)dt 2 + (1 - 2$)a 2 (t)5 ij dx i dx j , (25) 

where <3?(i, x) is the Bardeen potential which in this gauge coincides with the Newtonian potential. 
Defining the curvature perturbations IZ2 on comoving surface <j>2 = constant via 

n 2 = H-tl + $ , (26) 

and introducing the Sasaki-Mukhanov variables V2 

1 3 / 2 

v 2 = Z2K2 , z 2 = — j^- , (27) 

one obtains the following standard second order differential equation 

A + (k 2 c 2 s - Z f^V2 = Q. (28) 

Here the prime denotes the derivative with respect to the conformal time dr = dt/a{t) and c s is 
the sound speed associated with ^2 perturbations defined by 



Cs=l2 l = V 1 -/^- ( 29 ) 
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As mentioned before, we work in the ultra relativistic limit where 7/ 3> 1 correspondin g to a very 



small sound speed c s <C 1. This has interesting consequences for non-Gaussianities 



38, 



44]. 



To solve Eq. (|28p we need to find the form of Z2- Plugging different components of 22 from the 
background specified in the previous section, one obtains 

A ± + 2 

zf oc t A±+2 ocr A ±- x . (30) 

This in turn yields 

2 ^+ 5A ±+ 2 ^2 + ±)r-\ (31) 



z 2 ± (A| - 2A ± + 1) r 2 V A± 
where to get the final approximation the relation A± 3> 1 is used as explained in previous section 
Consequently, one can cast Eq. (1281) into more conventional form 



^ +(^-^-i^U T = 0, (32) 



where the indices /xj- are given by 

^4 + ^4 + 3 £± . (33) 

We have defined the "slow-roll" parameter e± = ^ ~ for the later convenience. As argued 
before, e± ~ 1/iVr < 1. 

One should solve Eq. (|32p for both t < t c and t > t c and glue the solutions via the matching 
conditions. In the limit where we can neglect the running of c s the solutions of Eq. (|32[) as usual are 
given in terms of Hankel functions H^{—c s kr) and H^{—c s kr). At the early stage of inflation 
when all physically relevant modes are inside the sound horizon, corresponding to —c s kr — > 00, the 



solutions are given in terms of positive frequency modes fc e lCakT created from vacuum. This 

r( 2 ), 



initial condition eliminates H^ l J.{—c s kT) and the solution of Eq. (|32h during the first inflationary 
stage is 

V2 _ = V^ e -(M-+l/2)/2^(D ( _ Csfcr) (34) 

After the collision, due to particles creation, the system is not in vacuum and both solutions 
H^2(~ c skr) and Hj^2(~ c skr) are permitted. We parametrize the solution after the collision in 



terms of the coefficients a and f3 via 

aH^(-c s kr) + (3H^(-c s kr)} . (35) 



V2+ = }Ell g ^(M+ + l/2)/2 



In the absence of any brane collision and particles creation a = l,/3 = 0, /x_ = //+ and the two 
solutions (134j) and (I35p coincides. Our job below is to find the coefficients a and j3 and read off 
the final value of curvature perturbations as a function of a and j3. For this purpose, we need to 
impose our matching prescriptions, joining Eqs. ([3l"j) and (f35|) at the time of collision r = r c . 

11 



B. Matching Conditions 



Imposing the matching condition is the most important step in our analysis of perturbations. 
There are standard prescriptions in literature for imposing the cosmological matching conditions 
in systems where there is a sudden jump in equation of state 



55, 



561 ] . More specifically, in models 



where the matching conditions are imposed on a comoving hyper-surface one requires that both 
the extrinsic and the intrinsic curvatures to be continuos across this hyper-surface. These result 
in the following matching conditions 



57 



59] 



m± = [n]± = o, 



(36) 



where <3? is the gauge invariant Bardeen potential and 1Z is the comoving curvature perturbations. 

In principle we can also impose these matching conditions in our system. However, as mentioned 
at the start of this section, we do not have a theoretical control on the dynamics of tachyon 
condensation, fields annihilation and particles creation. These shortcomings in turn result in a 
lack of knowledge of <3? and 1Z after the collision. To bypass these shortcomings we have chosen 
the phenomenological approach where only the perturbations of (j>2 are followed before and after 
the collision. In this approximation our model is effectively a single field scenario with the extra 
modification that the effects of particles creation are kept at the background level. This in turn 
results in different indices for Hankel functions, /x_ ^ In this view, our method is very similar to 



the method used in 



601 ] where the effects of phase transitions in a multiple field model is translated 



into a sudden violation of slow-roll parameters in an effective single field model. Therefore in our 
studies here, as in [60J, one expects that the matching conditions ([36]) are simplified into 



0, 



(37) 



i.e. both V2 and v' 2 are continuous across the surface of branes collision. 
Using the matching conditions (|3"7j) in Eqs. (fM|) and ([33]) one obtains 



a 



4 



H'W( x)H %(x)-H^(x)Hy(x 



r(2) 



A 2 ), 



(38) 



where the derivatives are with respect to the arguments of the Hankel functions and 



2^ 



(39) 
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IV. POWER SPECTRUM 



We are in a position to calculate the curvature power spectrum V-ji at the end of inflation. The 
curvature power spectrum is defined via 

Airk 

5 3 (k - k')Vn = — <ft(k')*^(k)) . (40) 

Since only the field 4>2 is present at the end of inflation we only need the power spectrum of 72-2 • 
Using Eqs. ([26]) and ([27]) one obtains 

O027 3 / 2 

At the later stage of inflation, when the mode of interest is outside of the sound horizon and 
—c s kr — > 0, one can approximate 

flg>(-q.*r) * -H^(-c s kr) ~ ^T( M+ ) (-^)^ + • (42) 

Plugging this into the power spectrum definition yields 

Vn = Vn W-a\ 2 , (43) 

where Vti represents the curvature power spectrum in the absence of branes collision 

Observationally Vn is fixed by the COBE normalization, Vn — 2 x 10~ 9 , when the mode cor- 
responding to the current Hubble radius leaves the horizon at about 60 e-folds before the end of 
inflation with k = aHj. One can use the COBE normalization to fix a combination of parameters. 
Specifically, using Eqs. (fT2|) and ([2Tj) we have 

Combined with Eq. ([2T|) the COBE normalization yields v^^l/^p — 2 x 10~ 3 . Assuming that 
A + < Nt = 60, this yiel ds A2 ~ 10 14 . This is the infamous fine-tuning problem associated with 
the DBI inflation 38, 6l}. Correspondingly m^jMp ~ 10 . 

To calculate the background curvature perturbations spectral index, n^, we note that at the 
time of horizon crossing dink = d In a + d In H + d In 7 = ( 1 — e+ + A f/^fH)H dt. Using this identity 
in ((SI) yields 
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a -PI 




FIG. 3: Here we plot the transfer function given by Eq. (j4"T)l with A_ = 45 and A + = 40. For large scales, 
x <C 1, it has a monotonic behavior as indicated by Eq. (|49p whereas for small scales, x 3> 1, it has an 
oscillatory behavior approaching the unity. The transition takes place for intermediate scales, x ~ 1, modes 
which leave the horizon around the time of branes annihilation. 



where to get the final result the relations H oc t 1 and <f>2 oc t 2 have 
ground. This indicates that is nearly scale invariant up to 0(e 2 , 



jeen used from the back- 
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43f]. 



The effects of branes collision and particles creation in the power spectrum Eq. ([430 are encoded 
in the transfer function \(3 — a\ 2 . In the absence of any feature, a = 1, (3 = so we obtain the 
standard single field DBI inflation 2^]. Below we would like to investigate the behavior of the 
transfer function \(3 — a\ 2 for different modes. Using the explicit formulae of a and j3 given in Eq. 
(f38[) one obtains 

rU/ \ 2 / \ 21 



1/9 |2 7F 



V J 



V Y 



x 



(47) 



As mentioned before we have defined x = —c s kr c = k/k c where k c represents the critical mode 
which leaves the sound horizon at the time of branes collision during inflation. One can check that 
in the absence of any feature, corresponding to fj,- = fi + , one obtains \/3 — a\ 2 = 1 as expected. 
As mentioned previously, the effects of branes annihilation and particles creations are such that 
Afj, ^ where Afi = fi + — 

The spectral index is also modified in the presence of the transfer function as follows 



nn - 1 



ax 



(48) 
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Now let us look at the shape of the transfer function for different length scales. In Fig. [3] we 
have plotted the transfer function. As a measure of the length scales of the perturbations, we note 
that x = 1 corresponds to the critical mode k = k c which leaves the sound horizon right at the 
time of branes annihilation. Then depending on whether the mode of interest leaves the sound 
horizon before the collision or after the collision one respectively has i <C 1 or j > 1. 

First consider the super-horizon scales, x <C 1. This corresponds to modes which leave the sound 
horizon before the collision. Using the small argument limit of the Bessel functions we obtain 

+ » > 2 (ir- <«> 

Our background is such that A + < A_ due to branes collision and energy transfer into radia- 
tion so A/x > 0. As a result the transfer function scales monotonically for super-horizon modes. 
Correspondingly, the change in spectral index is obtained to be 

An n = nTi-n\^ 2Ap . (50) 

We see from Eq. (j50|) that the spectral index for the super-horizon scales become more blue-tilted 
as compared to the background spectral index. This can also be seen in Fig. HI 

Now consider the sub-horizon modes, corresponding to x 3> 1. These are the modes which 
are deep inside the sound horizon at the time of branes collision. Physically one expects these 
modes to live in their local Minkowski background and should not be affected by branes collision 
and particles creation. Indeed, using the large argument limit of the Bessel functions J v {z) — 
^ cos (z-(v + ~)§ ) and Y u (z) ~ y^sin [z - (v + |)|) for z > 1 we find 

|/3 -a| 2 ~ l + 0(x- 2 ). (51) 

As expected the transfer function approaches unity so nn ~ n^. However, both the transfer func- 
tion and the spectral index are superimposed by oscillatory modulations with decaying amplitude 
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64] 



The most interesting effects of branes annihilation and particles creation is on the intermediate 
scales, x ~ 1. These are the modes which leave the sound horizon around the time of branes 
collision. In Figures [3] and [4] we have plotted the transfer function and the spectral index. As 
can be seen, for super-horizon scales the transfer function is monotonically increasing whereas for 
small scales it oscillates around unity. Correspondingly, the spectral index for super-horizon scales 
is well approximated by Eq. (I50p whereas for small scales it oscillates around n^. The non-trivial 
transition, happening at the intermediate scales, can be seen in both figures. 
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FIG. 4: Here we plot the spectral index n^. For super-horizon modes, x <C 1, our analytical formula Eq. 
(|48| gives a very accurate prediction for tin- For small scales, x 3> 1, the oscillatory behavior of is 
clearly seen as suggested by Eq. f|47[) where «r approaches its background value denoted by the dotted 
line. The transition takes place around x ~ 1 for modes which leave the horizon around the time of branes 
annihilation. The parameters are as in Fig. [3] 



As discussed before, the effects of branes annihilation and particles creation result in A/i 7^ 0. 
So far in the perturbations analysis we kept fi± and A± undetermined. Now let us specify Afj, in 
terms of our background model parameters. Starting with Eqs. JHJ), (f21~]) and ([33]) we obtain 

3 3 



2 



3 pi mf 
2A + p 2 m| 



(52) 

"2 

To get the final result we used the approximation that A + ~ A_. Physically, this is motivated 
from the assumption that the energy density associated with stack 1 is subdominant compared 
to the energy density stored in stack 2 so the energy transferred into radiation is much smaller 
than the background potential energy density. This assumption is necessary for the consistency of 
our matching condition. We have neglected the interference of perturbations 5<p\ and 5p r on the 
perturbations 5(f>2- As a consequence, we have imposed the simple matching condition (I37H . For 
this approximation to be valid, we require that p r /V{<f>2) < 1 or equivalently (A + — A_)/A + <C 1. 
We do not have good theoretical control on the masses mi and ni2 which were originated from 
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the back-reactions of Kahler moduli stabilization and background fluxes on the mobile branes. As a 
simple ansatz, we may take rrt\ ~ my- This corresponds to the case that the back-reactions imposed 
from the background compactification on a single brane is independent of its position inside the 
throat. This may not be an unreasonable assumption. With this approximation, the consistency 
of our approximation requires p\ <C P2- Also we are working in the probe brane approximations 
so the mobile stacks of branes should not deform the background AdS geometry significantly. For 
this to be the case we also need to impose p\ + P2 <C N. As an example, taking A + = 40, p\ = 1 
and P2 = 7 results in Ann ~ 0.01. 

One may also compare Ann induced from branes collision to the 0(e\) contributions to the 
background spectral index in Eq. ([15]) . We have A/x ~ 3(A_ — A + )/A+ ~ 3(A_ — A + )e^_. For a 
reasonable value of ( A_ — A + ) = 5 which was used in our numerical examples in Figures [3] and 0] 
we have A/i ~ 15e5_ or Ann ~ 30e^_. This indicates that the contributions of the branes collision 
to the spectral index is much bigger than the sub-leading 0{e\) contributions to the background 
n\. 

Q 

Observationally the result Ann > is disfavored [2] noting that the background spectral index 
is already nearly scale invariant as seen in Eq. (|46[) . To remedy this problem, it may be necessary to 
embed the simple model of ultra relativistic DBI inflation here to a generic setup of brane inflation 
where a period of slow-roll brane inflation takes place before DBI inflation starts. This can help 
to obtain a red spectral index for large scales whereas for smaller scales the spectral index can be 
blue as we have here. 



V. DISCUSSIONS AND CONCLUSIONS 



In this work we have studied the effects of branes annihilation and particles creation in models 
of DBI inflation where the branes are moving ultra relativistically. In a typical string theory 
compactification the process of branes and anti-branes annihilation is a generic phenomenon. It 
is therefore an interesting question how these collisions can affect the dynamics of inflation which 
may also be detectable in the sky. 

The process of branes annihilation is determined by the dynamics of open string tachyon con- 
densation. In a complicated background, such as the warped throat with background fluxes, this 
is a non-trivial process. The time-scale of tachyon condensation and branes annihilation is given 
by the inverse of the warped string mass scale, m s /i(ro) 1//4 . In our effective field theory approach 
the Hubble expansion rate is considerably smaller than this mass scale so for practical purposes 
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we can take the process of branes annihilation and particles creation to be instantaneous. In our 
set up the first inflationary stage is a two-field DBI inflation model. After the first stack of branes 
is annihilated with the background anti-branes a burst of closed strings modes are created. These 
particles are rapidly diluted by inflation and we eventually end up with a single field DBI inflation 
system. In order to connect the perturbations in the second inflationary stage to the perturbations 
generated before the collision we have to use proper matching conditions. Here are the major 
simplifications imposed in our analysis. It is assumed that the energy in field (pi, which is anni- 
hilated during inflation, is sub-dominant compared to energy stored in field cj>2 which carries the 
second stage of inflation. It is assumed that not only the background energy density in V\((j>i) 
is transferred into radiation p r but also the perturbations carried by 5<j>\ are also transferred into 
perturbations in radiation, 5p r . With these simplifications we can follow the perturbations in 02 
before and after the collision using the simple matching conditions (|37p . Although we did not take 
into account the perturbations 5<f>i and 5p r in our analysis but the effects of branes annihilation 
and particles creation are present at the level of background. In our analysis this translated into 
Afj, to be non-zero. Interestingly, this is very similar to model studied in 
ture is due to a sudden change in slow-roll parameters. In 
parameters was motivated from a phase transition 
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60, 



651 ] where the fea- 



651 ] the sudden change in s 



ow-roll 



661 ] during hybrid inflation model 
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69]. It 



is very interesting that our model, motivated from a completely different setup, shares the same 
results. 

We have calculated the curvature perturbations power spectrum. The effects of branes an- 
nihilation are encoded in the transfer function |/3 — a] 2 . For super- horizon scales, modes which 
are outside the sound horizon at the time of collision, the transfer function scales monotonically. 
Consequently, the spectral index is more blue compared to the background spectral index. On the 
other hand, for small scales, modes which are inside the sound horizon at the time of collision, 
the transfer function converges to unity with decaying superimposed oscillatory modulations. The 
non-trivial transitions in transfer function happen for modes which leave the sound horizon around 
the time of collision. The superimposed oscillatory modulations can have interesting observational 



consequences on CMB analysis 



65, 



711 ]. Furthermore, there are strong constraints from non- 



Gaussianity bounds on standard DBI inflation 
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431 ] . It would be interesting to calculate the level 



of non-Gaussianities produced in our setup to put further constraints on the model parameters. 
Our results here should be compared with the results obtained in 



was studied but in the context of slow-roll brane inflation. In 
the perturbations 5(j)\ are transferred into 5p r . However, in 



28| where similar problem 



281 ]. like here, it was assumed that 



281 ] the matching conditions (|36p have 
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been used instead of the simple matching conditions (|37p . As a result, the transfer function for 
arge and intermediate scales in [28] share the same behavior as in our case here. However, in 



281 ] (see also 70h72|] ) . the transfer function for small scales has a constant oscillatory modulation 



whereas in our case the amplitude of oscillations decays for these modes. For physical reasons, 
one expects that small scales should be blind to the process of branes annihilation and particles 
creation so the transfer function is expected to approach unity for these modes. To remedy this 
problem several methods were put forward in 28]. Our approach to use the matching conditions 



([3?]) may be interpreted as a pragmatic approach towards the methods speculated in [281 ]. 

To simplify our analysis we considered the model where the branes are moving ultra relativis- 
tically. In practice one expects that a successful brane inflation model consists both stages of 
slow-roll and fast-roll brane inflation. It is possible that during early stage of inflation the branes 
are moving slowly so one may get a few dozens of e-foldings during the slow-roll regime [17]. As 
the branes are moving towards the bottom of the throat they reach the speed limit and one should 
use the DBI inflation methods as we used in this work. The process of branes annihilation and 
particles creation can happen either during slow-roll regime, as in 28J, during the ultra relativistic 
fast-roll limit, as in our case here, or in between these two stages. In the third case, one can not 
perform analytical studies and a numerical investigation is necessary. It would be interesting to 
study the effects of branes annihilation in a generic model of brane inflation containing both stages 
of slow-roll inflation and DBI inflation. 
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